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We describe calculations of the properties of quantum fluids inside nanotubes of various 
sizes. Very small radius (R) pores confine the gases to a line, so that a one-dimensional 
(ID) approximation is applicable; the low temperature behavior of ID 4 He is discussed. 
Somewhat larger pores permit the particles to move off axis, resulting eventually in a 
transition to a cylindrical shell phase — a thin film near the tube wall; we explored this 
behavior for H2. At even larger R ~ 1 nm, both the shell phase and an axial phase are 
present. Results showing strong binding of cylindrical liquids 4 He and 3 He are discussed. 



1. Introduction 

The discovery of carbon nanotubes has provided a playground for theoretical physics 
analogous to that (~1970) based on the discovery of adsorption on flat, well- 
characterized surfaces. In the former case, excitement arises from the tantalizing 
possibility that one-dimensional (ID) physics can be tested by studying adsorbed 
gases near nanotubes, just as studies of monolayer films provided tests of 2D physics. 

Many groups have explored the properties of quantum fluids on the external sur- 
face of nanotube bundles and the interstitial regions within the bundles, stimulated 
by both the intriguing geometry and several experimental resultsP^EEEHHlour 
group has predicted several novel phase transitions for interstitial quantum fluids, 
including a high temperature (li quid- vapor) condensation and a BEC that exhibits 
4D (!) thermodynamic propertiesPElThis paper summarizes instead diverse results 
concerning quantum fluids inside single nanotubes, obtained with a variety of meth- 
ods. These studies are far from complete, with significant theoretical questions yet 
to be answered. The following section discusses the case of 4 He in ID, with appli- 
cations to small radius (i?) pores. Section 3 explores the behavior of absorbed H2 
as R increases, so that the ID approximation breaks down. Section 4 discusses the 
nature of films in large pores (R ~ 1 nm), where one encounters both a "cylindrical 
shell" phase of the film on the surface and a so-called "axial" phase, which is very 
much like the ID system in the small R case. 
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Throughout this paper we omit the details of both the adsorption potential and 
the techniques used in the calculations. Those can be found in existing publications, 
as well as a thesis and longer article currently being drafted™ Our emphases are 
new results, qualitative behavior, and significant open questions. 



2. Behavior in the ID Limit 

The ID 4 He system is interesting for several reasons. One is that the liquid is barely 
bound (by about 1.7 mK) and very rarefied (mean spacing about 2.7 ran!); in fact, 
the venerable Lennar d- Jo nes (LJ) pair potential is too weakly attractive to pro- 
duce this bound stateP^ A related fact is that the threshold interaction strength 
for ID binding of the liquid state coincides with that of the ID dimer!*" Recently, 
L.W. Bruch and C. Carraro (private communications) have shown that the cohesive 
energy of the ID many-body system has the ID dimcr binding energy as a lower 
limit, the two energies possibly coinciding. There is a closely related, intriguing 
aspect to the dimer problem. Consider the three-dimensional (3D) dimer problem, 
focusing on the s-wave channel. The radial Schrodinger equation for that problem 
coincides with the Schrodinger equation for the ID dimer. A key difference between 
D = 1 and D = 3 is the requirement that the wave function ip(r) vanish at 3D sep- 
aration r = 0. However, this difference is inconsequential for a ID system involving 
hard-core interactions. Hence, the 3D wavefunctions and spectra coincide, at least 
for the s-channel, with those of the ID problem. Putting all of this information 
together, it might be "expected" theoretically that all of the three energies agree, 
with the common value 1.7 mK. 

Recently, we hav e studi ed the thermal properties of ID liquid 4 He, using the 
path integral method.EIE] If one were to anticipate the behavior theoretically, one 
might treat the system with the Landau model, based on elementary excitations 
above the ground state. In the limit that the low-lying excitations are phonons, 
with T = speed s{p) at density p, this model predicts that the energy per particle 
AE (relative to the ground state energy) satisfies 

AE(T) = F(k B T) 2 , (1) 

F=^. (2) 
asp 

Preliminary results of the path integral calculations, in Fig.^ are consistent with 
this prediction at low T and high p; e.g., values of the coefficient F in Eq.Q]fit to the 
data at p = 2.5/nm (and higher) agree with the value predicted by Eq. |3 At p = 
2/nm and below, instead, the values of F begin to disagree and the departure from 
the quadratic dependence occurs at lower T. The latter is not surprising because 
the T 2 dependence of AE is expected only far below the ID Debye temperature, 
6 = hsp/(kBir)- Since O ~ 8 K at p = 2/nm and 9 ~ 30 K at p = 2.5/nm, different 
behavior of AE is expected for the two densities at 5 K. 

Departure from simple model behavior at very low p is also expected on general 
grounds, since the speed of sound becomes imaginary below the spinodal density, 
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Fig. 1. Energy per molecule of ID 4 He at p = 2.0 and 2.5/nm. Quadratic curves are based on 
the Landau-like model. (From Ref. HI) 

p ~ 0.245/nm (compared to an equilibrium density 0.36/nm). The phonon theory 
has no meaning at such low p. The observed departure from Landau-like predic- 
tions occurs at much higher density than that, however. This empirical behavior 
remains to be understood. At very low density, one might be inclined to use a virial 
expansion. Siber has evaluated the first virial correction to the ID 4 He system. 
His results show significant departure from the classical specific heat (1/2 Boltz- 
mann per atom) at relatively high T, even for very low density. This has nothing to 
do with exchange, which does not contribute because of the hard core repulsion.^ 
One expects higher order virial terms to contribute in addition, making this system 
particularly interesting to explore experimentally. 

3. Spreading Away From the Axis 

The T = properties of ID H2 were first calculated by Boronat, Gordillo, and 
Casulleras.^ Here, we report preliminary results of finite T behavior as the radius 
is increased from a very small value to R — 0.7 nm. 

Assuming a Lennard-Jones interaction between the gas adsorbate molecules and 
the inner pore surface, as adsorbate density is increased the gas will begin to adsorb 
onto the pore as a cylindrical shell film, at a distance from the surface on the order 
of the LJ parameter a. However, we expect that if the radius R of the pore is less 
than this characteristic distance, the gas-surface repulsion will heavily restrict the 
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transverse motion of the gas, and the two-dimensional (2D) shell will collapse into 
a ID line of molecules on the pore axis. 

Geometrically, this axial compression at small pore radii exhibits itself in the 
nanopore potential as a transition from a minimum near r ~ R — a and a maximum 
at r = (i.e., an off-axis potential well) when R > a, to a simple minimum at r = 
when R < a (an on-axis potential well). For the case of an infinitcsimally thin tube, 
this transition is analytically known to occur at R ~ 1.212<rP^We do not know an 
analogous analytic result for the case studied here, a pore within bulk material, but 
numerically the transition point is similar. 

For concreteness, we chose to study pores in MgO glass. The HVMgO LJ inter- 
action parameters are given by a — 2.014 A, e = 45.91 K .HI The H2-H2 interaction 
potential was taken to be of the Silvera-Goldman formal T he sys tem was stud- 
ied numerically via the path integral Monte Carlo algorithm! ^ I Temperatures 
between T = 0.5 K and 3 K, radii between R = 2 and 7 A, and densities of 
p = (2.55 x lQ- 3 /k)/R 2 and p = (7.13 x lQ- 2 /k)/R 2 . (The densities correspond to 
choosing a fixed number N = 1 or N = 28 of particles in a simulation cell of radius 
R and fixed length 125 A, with periodic longitudinal boundary conditions.) 

The results obtained were very similar for both densities and all ranges of tem- 
perature and radii studied. A representative plot of the H2 radial density distribution 
and the pore potential is given in Fig. |2Jfor radii between R = 2.5 and 3.25 A. 

The lowest pore size depicted, R = 2.5 A, is just beyond the transition point 
(2.45 A for the tubular case) to an off-axis minimum, where the H2 should switch 
from a ID axial phase to a 2D cylindrical shell. This is reflected in the potential 
energy curve (Fig. 0), which is very shallow for R = 2.5 A. Correspondingly, we see 
that at this radius the H2 molecules are still concentrated on the axis. However, with 
a remarkably small increase in pore radius, by 0.25-0.5 A, the distribution rapidly 
shifts to peak off axis, signifying the onset of the shell phase. It is interesting that 
the transition from axial to shell phases coincides in R with the shift of the classical 
potential minimum from on- to off-axis, even in the presence of quantum effects 
and interparticle interactions. One expects larger effects of interactions at higher 
density. 

It should be noted that the radii depicted in Fig. El are unphysically small for 
MgO pores of realistic size. We chose to study them despite this problem in order 
to theoretically study an axial to cylindrical shell transition in a pore geometry, 
which can only occur for small pores which are highly confining. In addition, there 
are some systems of larger radii in which the first layer to be adsorbed becomes a 
rigidly bound film, confining a fluid phase to a very localized vicinity of the pore 
axis Such possibilities help to justify the study of quasi-lD fluids. 



4. Large Pore Phenomena 

Relatively few simulation studies have been carried out for quantum fluids in "t 
ical" size nanotubes, R ~ 0.6 to 1 nm. Path integral calculations of Gatica et al\ 
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Fig. 2. H 2 radial probability density distribution (at T = 0.5 K, p = (2.55 X 10~ 3 /A)/_R 2 ) and 
MgO pore potential, for pore radii R = 2.5, 2.75, 3.0, and 3.25 A, as functions of dimensionlcss 
radius r/R. (Radial densities near r = are exaggerated due to finite size effects after normalizing 
the radial distribution by l/(27rr) to obtain the probability density.) 



reported the behavior of H2 over a limited range of R at T = 10 K. One of the more 
interesting phenomena is the pore-filling transition, shown for H2 in Fig. I^P^ At 
low /x, all of the molecules are localized within a thin layer, at r = 0.3 nm, located 
near the distance of closest approach to the nanotube. Above a threshold value of /x, 
the axial phase appears and grows rapidly with increasing fi. This axial phase can 
be thought of as an independent ID phase. A close analogy is the behavior of the 
second layer film of He or H2 on the surface of graphite, often treated by assuming 
that the only role of the first layer is to provide a holding potential!^ 

Analytical and numerical problems associated with matter in cylindrical geome- 
try are often computationally demanding, motivating the use of simplifying models 
that (we hope) capture the essential physics. Recently, we have explored such a 
description of the shell phase; the model assumes that all particles are constrained 
to lie on a cylindrical surface, r — R. One might expect that by varying R between 
R = and R = 00, one interpolates smoothly between ID and 2D behaviors. This 
is naive; instead, an intriguing "anomaly" arises: a significant enhancement of the 
binding occurs when the diameter of the cylinder, d — 2R, is comparable to the equi- 
librium separation r m j n in the pair potential. The condition d ~ r m m corresponds 
(for LJ interactions) to a/R ~ 1.7. Indeed, this argument does explain the R value 
corresponding to the maximum cohesive energy (seen in Fig. 0J of the "cylindrical 
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Fig. 3. Pore-filling transition of H2 in a tube of radius 0.6 nm, from Gatica et aili^l Results are 
densities as a function of r at 10 K and indicated chemical potentials. 



liquids" 4 He and 3 He. The 3 He case is perhaps the most dramatic, because its liquid 
state does not exist in either ID or 2D, while the cylindrical liquid 3 He is found to 
have cohesive energy as high as 1.26 K for R = 0.18 nm. 

These are variational results, obtained with Jastrow and Slater-Jastrow wave 
functions for 4 He and 3 He, respectively. Qualitatively similar, enhanced binding 
behavior was found for related problems involving similar binding problems on a 
cylindrical surface: He or H 2 dimers, a cry stalline lattice confined to a cylinder, and 
the virial coefficient of a classical fluidM^The origin of this general behavior is that 
two interacting particles can maximally exploit a divergent "specific area" when the 
interatomic separation is favorable. This occurs when the particles are on opposite 
sides of the cylinder, with separation |r 2 — ri| = d = r min . The specific area is defined 
as the cylindrical area residing within a separation interval [r, r + dr], divided by dr. 
For completeness, we note that analogous results for the dimer binding have been 
found by Aichinger et at, using both the simple mod el o f confinement on a cylinder 
and more realistic study of dimers inside a nanotubeU^ The optimal binding value 
found for R is very different in the two cases. 
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Fig. 4. Cohesive energy per atom of cylindrical liquids, as a function of R. Pluses are 4 He data 
(from Kostov et a/SHl) and circles are 3 He data (from C. Carraro, unpublished). 
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